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Abstract 

In this paper we present the results on Raman-active modes calculation in 
achiral {n,n) and (n, 0) carbon nanotubes (CNTs). vibrational spectra are 
derived as the eigenvalues of corresponding dynamical matrix at the F-point 
of Brillouin zone. Diagonal components and axy (= ctyx) are the only non- 
zero components of polarization tensor. Selection rules for Raman-active 
modes are determined by direct estimations of matrix elements responsible 
for the intensities of corresponding vibrational transitions. Few £^23 modes 
can be non-Raman-active because of their "fine" vibrational structure. We 
claim that numbers of Raman-active modes for achiral CNTs are: 5 ( (n, n) 
and (n,0), n-even); 6 ((n,n), n-odd); 8 ((n,0), n-odd). 

I. INTRODUCTION 

Carbon nanotubes (CNTs) are highly organized quasi-one-dimensional structures con- 
sisting of roUed-up graphene sheets. The variations in their diameters and helicities result 
in differences in electronic and vibrational properties fl]], [|2|, with Raman spectroscopy 
providing a powerful tool for their characterization. Currently, a few problems still exist in 
the calculation of the full list of Raman-active vibrations and their symmetry assignments. 
For some time it was believed that there were 15 or 16 Raman vibrations for achiral (n, 0) 
and (n, n) CNTs depending on whether n is odd or even 0]. The group-theoretical analysis 
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predicts 7 intensive Raman-active vibrations in the range up to 1600 cm~^ 0]. Decompo- 
sition of experimentally observed spectra into individual components was made following 
this predictions 0, 

Recently the total number of Raman-active vibrations was declared to be half of the 
previous estimates, namely 8 for both armchair and zigzag CNTs The corresponding 
calculations are based on the rod (or line) group theory worked out earlier for the 
stereoregular polymers P]. 

In the present work we perform a detailed analysis of the vibrational spectra of infi- 
nite CNTs with a view to address the unresolved issues. Our approach differs from the 
conventional approach, where the point and spatial symmetry groups are assigned first, 
followed by the full spectral analysis. Calculations are performed using Brenner's empirical 
potential for hydrocarbons [|l^]. Infinite CNTs are simulated by supercells with periodic 



boundary conditions. The full set of normal vibrational modes is obtained by the diagonal- 
ization of the dynamical matrix. The mode symmetries are assigned according to commonly 
accepted rules [|TT|, ||T^, [|T3|. Selection rules for Raman-active modes are found by an es- 
timation of matrix elements responsible for the intensity of the corresponding mode. The 
estimation is based on the direct comparison of the mode symmetry with the symmetry 
of polarization tensor components in polar coordinates. This approach does not require 
the preliminary complete group-theoretical analysis of CNTs. Nevertheless, the necessary 
symmetry assignment of vibrational modes is made. 

In this paper we present results for armchair (n, n) and zigzag (n, 0) CNTs with n = 
9 — 12. The detailed analysis of the vibrational modes demonstrates that there can be 
modes with identical symmetries that give different contributions to Raman spectra. It 
follows from the difference in the "fine" structure of these modes, and depends on whether 
vibrations of neighboring atoms are "in-phase" or "out-of-phase" . 

The manuscript is organized as follows: In the second section we present the details of 
the construction and analysis of the dynamical matrix. Sec. fTT] sketches the symmetry 
properties of achiral nanotubes. Sec. provides a brief description of the general rules 
for the estimation of the Raman-active vibrations. In Sec. we present results for achiral 



tubes (n, n) and (n, 0) where the vibrational spectrum of a (10, 10) nanotube is analyzed 
in more details as an example. Sec. |V| provides a summary of our results. 



II. CALCULATION OF VIBRATIONAL SPECTRUM USING THE 

DYNAMICAL MATRIX. 

The vibrational spectrum can be derived as the solution of the matrix equation: 

BX = Mu'^BX, (1) 

where D is the dynamical matrix and B is the matrix of coefficients of quadratic form of 
kinetic energy Ekm = ^ bij Vi vj, {vi, vj are generalized velocities of atoms i and j). {uk} 
and {Xk} are, respectively, the eigenvalues and eigenvectors fully describing the vibrational 
spectrum. M is the mass of carbon atom. Matrix elements Dj^ are given by: 

where E is the potential energy of the system; ctj, j3j are generalized coordinates of atoms i 
and j. In Cartesian coordinates B in (|l]) is the unit matrix and vibrational spectrum is ob- 
tained by diagonalization of matrix D. Corresponding eigenvectors contain all information 
on symmetry and relative amplitudes of atomic displacements from the equilibrium. The 
Brenner's potential [|1^] is used for modeling nanotube structure and vibrational properties. 
Analytical form of this potential allows to get second derivatives analytically. Cylindrical 
coordinates greatly facilitate the construction of the dynamical matrix: it is sufficient to 
calculate the matrix elements for two neighboring atoms and then the dynamical matrix 
is reconstructed by symmetry. The aspect ratio of CNTs is very large and it allows to 
consider the CNTs as infinite. We use a supercell with periodic boundary conditions to 
model the infinite tube. 

If a supercell consists of atoms, the rank of the dynamical matrix is 3 x 3 and 
Eqn. d^) has 3A^ solutions giving the full set of normal vibrations. We consider modes 
only in the center of the Brillouin zone (F-point with /c^ = 0). Modes with /c^ 7^ give 



contribution to Raman-spectra for finite CNTs |14], [15], and will be ommited in the present 



consideration. As a result, the rank of the dynamical matrix decreases to 3 A^' x 3 A^', where 
3 A^' is the number of vibrations with k-. = 0. 



The symmetry assignments are standard (see, e.g., ||rT|, [|T2|, [|T^). For one-dimensional 
representations letters A and B stand for symmetric and antisymmetric modes with respect 
to symmetry operations; subscripts '1' and '2' - denote symmetric and antisymmetric 
modes with respect to the rotation about two-fold axes perpendicular to the principal axis. 
The letters ^g' and 'm', - symmetric and antisymmetric modes with respect to the inversion. 
In two-dimensional representations Ek, the subscript k stands for modes which transform 
to itself by the symmetry operation. An advantage for the choice of S2n as the principal 
axis over the C„ axis for CNTs is given in Sec. 0. The results are summarized in Tables I, 
II and will be discussed below. 

There are A and B modes with "double" symmetries with respect to the rotation about 
C2 axes perpendicular to the principal axis S2n- An example of such (Raman- inactive) 
mode in {n,n) CNT with even n is given in Fig. 1. The mode shown is the Bg tangential 
mode where all atoms in every layer0 have equal displacements. There are three classes of 
C2 rotational axes perpendicular to the tube axis (C2, C2, C2) for even n. (More detailed 
analysis of CNT symmetry elements will be given in the next Section). One can see that 
mode Bg is antisymmetric with respect to €'2 and €'{ axes and is symmetric with respect to 
C'2 axis. This is the reason why subscripts '1' or '2' are omitted for some A and B modes 
in Table I. 

Following these rules, the symmetry analysis shows that the vibrational modes of (n, n) 
and (n, 0) CNTs with even and odd n have following symmetries: 



All atoms with equal z-coordinates form a layer. 
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= 2 Ai, + 2 + 2 A„ + 1 5i„ + 1 + 4 Bg+ 

n—l)u 

^VS'' = 2 + 2^23 + 2^3 + + 152„ + 45„ + 6Ei„ + 6^25 + . . . + 6E(„_i)3 

rn-e.en = 2 Ai, + 1 ^2, + 3 A„ + 1 5i„ + 2 + 3 Bg+ 

3 -Elg + 3 Elu + 3 £'29 + 3 E2U + . . . + 3 E{^n-l)g + 3 -E'{n-l)n 

pn-odd _ 2 Ai, + 1 ^2, + 3 + + 252„ + 3fi„ + 6Ei„ + 6^23 + . . . + 6E(„_i), 

(3) 

Four zero-frequency modes are not shown in @. The rotation about z-axis has A2g sym- 
metry (see Fig. 2a); the translation along the 2;-axis has i?2u symmetry i; translations in 
(x, y) plane are doubly-degenerated and have Ei^ symmetry (Fig. 2c). Symmetries of zero- 
frequency modes differ from symmetries given in because we use S2n as the principal 
axis. 

Before proceeding with the analysis of the Raman-activity of calculated vibrational 
modes we briefly recall the main symmetry elements in CNTs. 

III. SYMMETRY ELEMENTS IN CNTS. 

An analysis of the vibrational modes and comparison with polarization tensor com- 
ponents requires a knowledge of necessary symmetry elements of CNTs. Note that the 
symmetry analysis is independent on the form of interatomic potential; this choice deter- 
mines only frequencies but not symmetries. 

The full group-theoretical analysis of CNTs is beyond the scope of this work i. A list of 



symmetry is assigned because this mode is antisymmetric with respect to 52„ operation which 
is a rotation by an angle 27r/2n and reflection in the horizontal plane passing through centers of 
the slanted bonds. If the C„ axis is taken as the principal axis then this mode has A2u symmetry 
since it is symmetric with respect to the rotation (Fig. 2b) 

^One can find the group-theoretical analysis in Q, and nonsimmorphic rod-group is: 
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allowed point symmetry elements for (n, n) CNTs with even and odd n is given in Table III. 
Symmetry elements for Dnh and Dnd point groups are also shown for comparison. 

Fig. 3 illustrates symmetry operations for (n, n) CNTs with even and odd n (n = 6 in 
Fig. 3a, and n — 5 in Fig. 3b). Two different forms of unit cells are used, since different 
symmetry elements are clarified in different representations. Two top panels show chains of 
hexagons chosen as unit cells. Top panel in column (b) (odd n) shows a two-fold rotation 
axes, C2, passing through centers of a hexagon and an opposite horizontal bond. This 
axis splits into two classes of two-fold rotation axes for even n, as shown in the top panel 
of column (a). One rotation axis, Cg, passes through centers of opposite hexagons, while 
the other, , passes through centers of opposite horizontal bonds. The (n, n) CNT with 
even n is thus associated with two classes of reflection planes: a'^ and a" passing through 
C'2 and C2 axes, respectively. By contrast, there is only one class of reflection planes, a[, 
associated with a CNT with odd n. Moreover, while the point i in the top panel for even 
n is an inversion center, the corresponding point O' in the right panel is not an inversion 
center. Two top panels also show that both tubes have horizontal reflection planes ah 
passing through centers of hexagons. 

In order to illustrate other symmetry elements we show the armchair fragments as the 
unit cell for both tubes in middle panels in Fig. 3. One can see that there are additional 
two- fold rotation axes, Cg", passing through centers of opposite slanted bonds. Additionally, 
both tubes have 2 n-fold roto-reflection axes, 5'2„. The inversion center in the case of odd 
n (right panel) is i and inversion is absent for even n. 

The bottom panels in Fig. 3 are the side view of both tubes with corresponding sym- 
metry elements. 

To summarize: there are symmetry elements, C„, S2n, Cv, C2"; well as two types of 
fixed points, i (inversion) and O, which are common to CNTs with both even and odd n. 
Symmetry elements C2, C2, cr^, cr'^ are inherent to CNTs with even n, while elements C2 



Cnv X S2n ■ 
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and a[ are characteristic to odd n. 

This analysis also demonstrates that all symmetry elements of both Dnd and D^h groups 
are permitted in CNTs (see Table III and Fig. 3). One can say that CNTs 'combine' 
symmetry properties of Dnh and Dnd point groups. In particular, CNT has three classes of 
two- fold axes for even n: (Cg, Cg, C2). These axes are the combination of two classes for 
Dnh symmetry group {C'2, axes) and one class of C2 axes for Dnd {C2 axes). 

In the case of zigzag CNTs an analysis is similar. Thus, with a knowledge of CNTs 
symmetry elements and the calculated vibrational spectra we can identify Raman-active 
modes in CNTs. 

IV. HOW TO REVEAL RAMAN- ACTIVE MODES? 

In this Section we recall main rules for determining Raman-active modes. Let the 
transition between states Ei and Ej is associated with wave functions ipi and ipj. Then the 
probability for the ipi ipj process is proportional to the square of the transition moment 
Mij, where 

Mij = I r^M^i dV. (4) 

In Raman process, M is the dipole moment induced by the external field E. The induced 
dipole moment is P = aE, where a is a polarization tensor. Thus, the selection rules in 
Raman spectra are determined by the values of integrals 

j r.a^^'i^^dV, (5) 

where a^^i is a component (or a linear combination of components) of the polarization 
tensor (,^, ^' = x,y,z). The transition i ^ j is allowed if (|) is nonzero for at least one 
component a^^/. The integration is performed over the full set of electronic and nuclear 
coordinates. Direct calculation of (^ is very complicated problem. So, the usual practice 
is to consider the electronic, vibrational and rotational degrees of freedom separately. Thus 
the wave function can be written as the product ip = ip^^ ip^^^ tp'''"^ . As a result only the 
following integrals should be estimated for Raman transitions: 
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j^fa^^.^fdV. (6) 

vibrational wave function ip^^^ is real and is represented by an eigenvector of the dy- 
namical matrix. Since we consider the F-point of Brillouin zone, where the dependence on 
the z coordinate is absent, the space integration in (^) is replaced by the area integral 

j ^fa^^,^T'dxdy. (7) 

It is convenient to return back to the polar coordinates (^, ^' = r, (p). In the simplest 
case the ground-state vibrational mode ipi is fully symmetric and is a constant in polar 
coordinates. Moreover, the vibrational wave function is the discrete function (eigenvector 
of the dynamical matrix) . Then (0) transforms to the sum: 

2n 

E v^r^e^', (8) 

where the summation is over atoms m, m = 1, 2, . . . 2 ri of only one layer, and products of 
polarization tensor components with m-th elements ipY' of eigenvector of j-th normal mode 
are calculated for angles 0*", corresponding to polar coordinate of atom m. This approach 
is very convenient as both polarization tensor components and eigenvector elements of 
vibrational modes have very simple representation in polar coordinates for CNTs. 

As is demonstrated in Appendix, only a^x, c^yy, c^xx ± Ciyy and axy components of 
polarization tensor differ from zero for CNTs and are shown in Fig. 4. 

Many one- dimensional modes in can be excluded from candidates for Raman vi- 
brations, (i) All modes with subscript 'u' are eliminated as all components of polarization 
tensor are symmetric with respect to this operation; (ii) All B modes are excluded as an- 
tisymmetric with respect to operation, (iii) modes are also non- Raman- active as 
they are asymmetric with respect to C2 operations. Hence, only two fully symmetrical Aig 
modes are Raman-active. 

Doubly-degenerated E^g modes in polar coordinates behave as sin(/c0 -|- ip), where ip - 
some phase. Their schematic view is given in Fig. 5. The signs assignment is arbitrary. 
For instance, the clockwise direction for tangential displacements and radial displacements 
out of center have sign. 
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Comparison of Fig. 4 and Fig. 5 results in the fact that only mode can be Raman- 
active. Really, the shapes of components of polarization tensor coincide only with the shape 
of mode, and signs of their lobes alternate in the same manner. Actually, it is necessary 
to estimate expression (|^), which can be reduced to 

/ sm{2 (j) + e) sm{k (j) + d(j), (9) 

^0 

and is different from zero only if = 2. (sin(2 + 6) is the general representation for axy 
and axx-oyy components of polarization tensor, and phases ip - for vibrational modes and 
6 - for polarization tensor are related as their principal axes coincide.). It means that only 
modes can be Raman-active. The conclusion that only Aig and E2g modes are Raman- 
active is valid for all types of achiral CNTs, both armchair and zigzag with arbitrary n (see 
Table II). 

V. RAMAN-ACTIVE MODES IN ACHIRAL CNTS 

Two Aig modes for all types of CNTs are Raman-active as they are fully symmetrical 
"breathing" and tangential modes (see Fig. 6). Now we will demonstrate that one of E2g 
modes is not Raman-active on the example of (10, 10) CNT. All E2g modes are shown in 
Fig. 7. One can see the difference between mode (c) and three other modes, this deference 
being in 'signs' of relative displacements of neighboring atoms. These displacements are 
predominantly "in-phase" for (a), (b) and (d) modes and are "out-of-phase" for (c) mode. 
Hence, these "phase" relations of different modes give different contributions to (P): this 
sum is not equal to zero for modes (a), (b), (c) and is zero (or very close to zero) for mode 
(c). It means that mode (c) is not Raman-active, and Raman spectra of (10, 10) CNT have 
5 Raman-active vibrations (see Table II). These conclusion is valid for all (n, n) CNTs with 
even n. 

Situation with E2g vibrations for (n, n) CNTs with odd n is analogous. Among six E2g 
modes there are two "out-of-phase" modes which are not Raman-active, and total number 
of Raman-active vibrations in this type of CNTs is 6. 
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There are no "out-of-phase" modes in (n, 0) CNTs, and the total number of Raman- 
active vibrations is the sum of Aig and mode types: 5 Raman-active modes for (n, 0) 
CNTs if n is odd, and 8 Raman-active modes if n is even. 

The reason why armchair CNTs differ from zigzag CNTs consists in the fact that the 
latter have one atom in the primitive cell, and the former - two atoms. Two atoms in the 
primitive cell for armchair CNTs allow intrinsic vibrational "fine" structure. Qualitatively 
it can be demonstrated as follows. 

We consider Ekg modes, and the armchair ring is chosen as the unit cell. This unit 
cell consists of 4 atoms from two layers (see Fig. 8a, where (6, 6) CNT is shown, as an 
example). Let us open this fragment into a planar structure (Fig. 8b) and further project 
into a linear chain (Fig. 8c). There are two atoms in the primitive cell (two adjacent 
atoms in a layer) since axis was chosen as a principal axis. So every filled circle in 
Fig. 8c represents one primitive cell of two atoms. This gives a chain with cyclic boundary 
conditions consisting of ric = 2 n units, and its solution is convenient to write in the form: 

u{k, I) = a cos{2tt k I /ric), (10) 

where u{k, I) is a deviation of the unit / for k-th type of normal vibration (A; = 0, 1, . . . ric—l), 
and / numerates the units in the chain (/ = 1, 2, . . . ric)- 

So, the full vibrational spectrum of (n, n) CNT in the center of the Brillouin zone can be 
approximated by a solution of the one-dimensional chain with cyclic boundary conditions 
and some effective potential. 

The solution ([T0| ) gives all one- and two-dimensional modes. Really, if k = then all 
units have equal deviations u{0, I) = a of the same sign, which corresponds to A modes. If 
k = nc/2 then all units have equal absolute values of deviations, but the signs of relative 
deviations alternate u{nc/2, I) = (— l)'a, which corresponds to B modes. 

Modes Ekg are formed by pairs of solutions (p!0|) with 'conjugated' values of k: k and 
He — k. Really, 

a cos(2 TT {ric — k)l/ Uc) = a cos(2 it kl/ ric) , 
and it gives (n — 1) E^g modes for (n, n) CNTs in agreement with (||). 
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The solution ([T0|) for = 2 is shown in Fig. 8c and corresponds to mode of the 
nanotube. This mode can be represented schematically in polar coordinates as shown in 
Fig. 8d and is equivalent to mode E2g in Fig. 5. 



VI. CONCLUSIONS 

The main goal of the present communication is the calculation of numbers of Raman- 
active modes in achiral CNTs {n, n) and {n, 0) {n = 9-12). vibrational modes are calculated 
in F-point of Brillouin zone, and Raman-active modes are identified by an estimation of 
matrix elements responsible for the intensity of corresponding transitions, and no prelim- 
inary full group-theoretical analysis of CNTs is necessary in this case. For achiral CNTs 
only three diagonal and a^y = <yyx components of polarization tensor are not equal to zero. 

We have shown that there are one or two modes for (n, n) CNTs which are not 
Raman-active. It follows from the fact that the primitive cell of these CNTs has two 
atoms, and vibrational modes have "in-phase" and "out-of-phase" types of vibrations, the 
latter being non- Raman- active. The primitive cell of {n, 0) CNTs consists of one atom and 
these CNTs has no analogous complications. 

The calculated numbers of Raman-active modes are five or six for {n, n) CNTs for even 
and odd n, correspondingly. For (ra, 0) CNTs there are five Raman-active modes for even 
n and eight Raman- active modes for odd n. 

Our primary goal was to formulate the selection rules for Raman-active modes and 
less attention was paid to calculations of mode frequencies. More precise values of mode 
energies can be obtained using more sophisticated potentials, it e.g. Menon's tight-binding 



scheme |jT6|, [O]. 
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APPENDIX A: POLARIZATION TENSOR COMPONENTS FOR CNTS 



1) We consider non- resonance Raman transitions. For simplicity achiral CNT are rep- 
resented as homogeneous hollow cylinder of sufficiently large radius. 

2) Let external electric field is applied along the x-axis perpendicular to the vertical 
nanotube axis z. This field can not induce polarization along 2;-axis as the "top" and the 
"bottom" of CNT are equivalent. By symmetry the same is valid for |/-direction of applied 
field. Hence, a^z = cXzx = Ciyz = Cizy = 0. 

3) Now we define the angular dependence of components a^x and a^y in polar coordi- 
nates. External electric field is aligned along the x-axis and we consider the {xy) section 
of CNT (see Fig. 9) . Polarization vector at point A, (0 = 0) has only normal component: 
Pn = otnE , and a„ is the normal polarizability. There is only tangential component at 
point B {(f) = tt/2) and Pr = (Xr E, where a,- - tangential polarizability. Obviously, a,- > an 
as polarizability in the plain is larger compared to perpendicular direction. 

4) Consider now point C at an angle (p to a;-axis. Electric field at this point has 
normal {En = E cos 0) and tangential {Er = E sin 0) components. Normal and tangential 
components of polarization vector at this point are: 




(Al) 



and X- and ^/-components of polarization vector P' at point C are: 




(A2) 



Substituting in system (^) values for P^ and from ( [Al| ) one gets: 



P^ = {an cos^ 4> + ar sin^ 0) E 

Py = (ttn sin (p cos (j) — ar sin (p cos (p) E 



(A3) 



And it follows that 



a 



XX 



an cos^ (p + ar sin^ (p 



(A4) 



< 



a 



xy — 



i {ar — an) sin 2 (p 
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One can get expressions for other non-zero components of polarization tensor in polar 
coordinates: 



aj- cos^ + q;„ sin^ 

ar + an (A5) 

(q!„ — a-r) COS 2 (f) 

Schematic view of all non-zero components of polarization tensor is shown in Fig. 4. 



(Xyy 

O-xx ^yy 



13 



REFERENCES 

[1] M.S.Dresselhaus, G.Dresselhaus, and P.C.Eklund, Science of FuUerenes and Carbon 
Nanotubes (Academic Press, San Diego, 1996). 

[2] R.Saito, G.Dresselhaus, and M.S.Dresselhaus, Physical Properties of Carbon Nanotubes 
(Imperial College Press, London, 1998). 

[3] R.A.Jishi, L.Venkataraman, M.S.Dresselhaus, and G.Dresselhaus, Chem.Phys.Lett., 
209, 77 (1993). 

[4] R.Saito, T.Takeya, T.Kimura, G.Dresselhaus, and M.S.Dresselhaus, Phys.Rev.B, 57, 
4145 (1998). 

[5] A.M.Rao ct al. Science, 275, 187 (1997) 

[6] C.Journet et al.. Nature (London), 388, 756 (1997). 

[7] O.E.Alon, Phys.Rev.B, 63, 201403R (2001). 

[8] M.Damnjanovic, I.Milosevic, T.Vukovic, and R.Sredanovic, Phys.Rev.B, 60, 2728 (1999). 
[9] I.Milosevic, R.Zivanovic, and M.Damnjanovic, Polymer, 38, 4445 (1997). 
[10] D.W.Brenner, Phys.Rev.B, 42, 9458 (1990). 

[11] R.L. Flurry Symmetry groups. Theory and chemical applications, (Prentince-Hall, INC, 
New Jersey 1980). 

[12] M.M.El'yashevich Atomic and molecular spectroscopy, (Editorial URCC, Moscow, 2001, 
in Russian). 

[13] E.B.Wilson, J.C.Decius, P.C.Cross The symmetry of infrared and Raman vibrational 
spectra, (McGraw-Hill Publishing Company LTD, New York, London, Toronto 1960. 

[14] R.Saito, T.Takeya, T.Kimura, G.Dresselhaus, and M.S.Dresselhaus, Phys.Rev.B, 59, 
2388 (1999). 

[15] E.Richter, K.R.Subbaswamy Phys. Rev. Lett., 79, 2738 (1997). 

[16] M.Menon, K.R.Subbaswamy, and M.Sawtarie, Phys.Rev.B, 48, 8398 (1993). 

[17] M.Menon, E.Richter, and K.R.Subbaswamy J.Chem.Phys., 104, 5875 (1996). 



14 



TABLES 

TABLE I. The vibrational spectrum of (10, 10) CNT calculated using the Brenner potential. 

The mode frequency lo (in cm~^) and its symmetry S are given in the first and second columns. 
The polarization tensor elements a^^' are also shown. Raman-active vibrations are underlined. 
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TABLE II. Raman-active modes of (n, n) and (n, 0) CNT for n = 9 — 12 calculated using the 
Brenner's potential. 
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JTl cqUcIiClco lOIli J cLilLl oyiliilicUlico Ui XLdlllcliIl ctCLlVc IliULlco 




J^O \^2g)i \-t^lg)-i ODD {-t!j2g), IDDO \i^2g)-i IDD^ l,-C/25j) iDe4 l^^lgj 


(^(^ ^^^^ 








(12,12) 


9 (E20), 125 (y4i„), 278 (EaJ, 1680 (£;2o), 1689 (^i.) 


(9,0) 


47 iE2g), 286 (A13), 362 (^23), 591 {E2g), 741 (£;25), 1144 (^29), 1672 (^29), 1677 (^i^) 


(10,0) 


38 iE2g), 258 (A15), 542 (E25), 1676 (£;29), 1680 (^i^) 


(11,0) 


31 {E2g), 235 (^15), 348 (^25), 499 (^25), 695 {E2g), 1147 (^2^), 1679 (£;2p), 1683 (^i^) 


(12,0) 


26 {E2g), 216 (Aig), 462 (^2^), 1681 (£;2<;), 1685 (A^g) 
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TABLE III. Symmetry operations for Dnh and D^d point groups and symmetry operation for 
{n,n) CNTs. Here E is the unit operation; is the rotation about the z-axis by angle 27rm/n; 
6*2, 6*2, C2' are rotations about (72,6*2, C2" axes (see Fig. 3) by an angle n; ah is the reflection 
in the plane perpendicular to the main symmetry axis; (t„ is reflection in the plane including the 
main symmetry and one of C2 axes; is reflection in the plane which includes the principal axis 
and bisects the angle between two two-fold axes; 5™ is the rotation about the principal axis by 
angle 2iTm/n and reflection in the plane ah'- S^^ = cr/j; i - is the inversion operation. 













J^nd 




(n, n) 


(n,ri) 








even n 


odd n 


even n 


odd n 


even tt, 


odd n 








1 


1 


1 


1 


1 


1 


n ' 


= 1,- 


. n — 1 


+ 


+ 


+ 


+ 


+ 


+ 




= 1,.. 


.n-1 


+ 


+ 






+ 


+ 




= 1,3. 


.2n-l 






+ 


+ 


+ 


+ 




ay 




+ 


+ 






+ 


+ 








+ 


+ 






+ 


+ 








+ 




+ 


+ 


+ 


+ 




C'2 








+ 


+ 


+ 


+ 








+ 


+ 






+ 










+ 








+ 


+ 




i 




+ 






+ 


+ 


+ 
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FIGURE CAPTIONS 



Fig. 1. An example of Bg mode (a; = 1172 cm~^) with "double" symmetry is shown as a 
side view of an armchair CNT with even n. The mode is antisymmetric with respect to 
rotation about C2 and C2 axes and is symmetric with respect to rotation about C2 axis. 
Axes are perpendicular to the principal axis pass through centers 

of opposite horizontal bonds, centers of opposite hexagons, and centers of opposite slanted 
bonds, respectively. 

Fig. 2. Zero- frequency modes: (a) - rotation about the principal axis; (b) - 
translation along the principal axis; (c) Ei^ - translations along x and y axes perpendicular 
to the principal axis z. 

Fig. 3. Unit cells of (6, 6) (left column) and (5, 5) (right column) CNTs with their symmetry 
elements. Top panels show unit cell chosen as a chain of hexagons; middle panels present 
armchair unit cells. The side view of tube fragments is given in the bottom panels. 

Fig. 4. Schematic representation of non-zero polarization matrix elements in (x, y) section of 
CNTs: — cin cos^ 0-1-0^ sin^ 0, ayy — a-r cos^(l)+an sin^ 0, a^y — ocyx — \ (a„— a,-) sin 20, 
0!.xx — OLyy — {an — Qr) COS 20. Note the sign alternation for axy and axx — ocyy components. 

Fig. 5. Schematic representations of Eig - Egg modes. Note the sign alternations in lobes 
of vibrational modes, representing the relative displacements of atoms from equilibrium 
positions. 

Fig. 6. Top view (in the (x, y)-section of (10, 10) CNT) of Raman-active Aig modes: (a) 
"breathing" mode (a; = 150 cm~^); (b) high-frequency tangential mode {u — 1686 cm~^). 

Fig. 7. All E2g modes for (10, 10) CNT in the {x,y) plane: (a) uu = 12 cm"^; (b) uu = 331 
cm~^; (c) u! = 826 cm~^; (b) uj = 1673 cm~^. 

Fig. 8. (a) Armchair unit cell of (6, 6) CNT; (b) unit cell projection on the plain; (c) 
Schematic representation of the unit cell as a linear chain (with cyclic boundary conditions) , 
where every pair of atoms of one layer is substituted by one unit (the vibrational mode for 
A; = 2 is shown); (d) representation of this mode in polar coordinates. 

Fig. 9. Components of electric field and polarization vectors in the cross-section of CNT 
for different polar angles (p. 
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